Manipulation of double-dot spin qubit by continuous noisy measurement 
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We consider evolution of a double quantum dot (DQD) two-electron spin qubit which is continu- 
ously weakly measured with a linear charge detector (quantum point contact). Since the interaction 
between the spins of two electrons depends on their charge state, the charge measurement affects 
the state of two spins, and induces non-trivial spin dynamics. We consider the regimes of strong 
and weak coupling to the detector, and investigate the measurement-induced spin dynamics both 
analytically and numerically. We observe emergence of the negative-result evolution and the system 
stabilization due to an analog of quantum Zeno effect. Moreover, unitary evolution between the 
triplet and a singlet state is induced by the negative-result measurement. We demonstrate that 
these effects exist for both strong and weak coupling between the detector and the DQD system. 



I. INTRODUCTION 



Quantum dots (QDs) recently have attracted much at- 
tention as very suitable candidates for studying evolution 
of a single quantum system. Moreover, QDs are very 
promising candidates for quantum information process- 
ing. An electron spin in a quantum dot is a good repre- 
sentation for a qubit, being a natural two-state quantum 
system. Manipulations of the single QDs and double-QD 
(DQD) systems can be achieved using external magnetic 
fields and gate voltages^, which implement injection of 
an electron (or pair of electrons for DQD) in a given 
spin state, or some unitary rotation in a two-spin space. 
However, it is not easy to achieve the full set of transfor- 
mations in the Hilbert space of two coupled spins^^: 
some transformations require a strong gradient of mag- 
netic fields on a nanometer scale, and advanced tech- 
niques are needed for performing spin rotations rapidly 
and reliably. 



Thus, it is interesting to explore other possible ways 
of manipulating electron spins in DQD structures. In 
particular, it is natural to study whether measurement 
of the charge state of the DQD may help drive the de- 
sired evolution of the two-spin system. The weak contin- 
uous measurement^^ 9 -, which monitors the system in 
question, and therefore affects its evolution, may provide 
an additional and useful tool for controlling the electron 
spins in DQD systems. At some level of description of 
a probed quantum system plus detector, it is postulated 
that the fundamental measurement process consists of 
direct particle detections (e.g., photons, electrons, etc.) 
or absence of any detection. The absence of a detection 
for specific time interval constitutes a negative result. In 
both cases ("detection" or "no detection") the quantum 
state of the system evolves according to the provided in- 
formation and the evolution can be described using quan- 
tum Bayesian inference. 9 E.g., for a two-level atom which 
did not decay by time t the evolution of the atom's den- 
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FIG. 1: Schematic of a DQD spin system measured by a 
QPC; the QPC current I(t) is affected by the charge state 
of the DQD. In quantum Zeno regime the tunneling between 
dots is suppressed while the spin qubit undergoes continuous 
negative-result evolution. 
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where e _rapt is the conditional probability not to decay 
from the excited state |e) (ground state \g) is supposed 
stable), T sp is the spontaneous decay rate, and J\f = 
p gg (0) + Pee(0) e _rspt . The classical Bayesian update of 
the diagonal elements, Eq.([T]), reflects an information- 
related evolution. Also, Eq.Q shows that the coherence 
ratio | Peg | / 'v/ 'Pee Pgg is conserved; in particular, a pure 
state remains pure. 

The above evolution has been proposed and discussed 
mainly in the context of quantum opticsii^'ii^ For solid 
state systems, it has been realized recently, using an arti- 
ficial two-level system: a superconducting "phase" qubit, 
measured via tunneling of its superconducting phase. 11 
The detection of tunneling was considered as an instan- 



taneous process 



in the highly non-linear detector 



(dc-SQUID) that implies detector relaxation rates much 
faster than the qubit "decay rate". In this paper we 
consider another realization of a solid state qubit and a 
linear detector measuring continuously. Nevertheless, it 
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can behave similarly to the above examples in certain 
measurement regimes. We will consider a semiconductor 
double quantum dot (DQD) spin system (FigQ]) where 
two electrons can be in either of the dots and the rele- 
vant states are spin-charge states ^ 13 ' 14 ! 15 At relatively 
small time scales (from ns to /is) the electron tunneling 
between the dots preserves the system's total spin that 
allows spin-selective evolution and leads to the so called 
spin-to-charge conversion^ This allows the measurement 
of the DQD spin system by a charge sensitive detector 
such as a quantum point contact (QPC ) 15 ' 16 or single 
electron transistor (SET)f^ 

For our system (FigfTJ) the spin qubit 15 is formed by the 
singlet spin state |5(1, 1)) and the triplet state |T (1, 1)) 
that has a zero spin projection on an external magnetic 
field (here, e.g., (1,1) denotes a charge state when one 
electron is in the left dot and the other is in the right dot). 
Spin conservation allows transition between the singlet 
|5(1, 1)) and the localized singlet, 15(0, 2)), while a tran- 
sition from the triplet state is spin-blocked. Thus, the 
state |5(0, 2)) can be used for preparation and subsequent 
measurement of the spin qubit. For a weak measurement 
the system-detector coupling V is finite and defines the 
time scale of measurement evolution (in particular, 1/r 
is of the order of the time needed to distinguish the two 
different charge states by the QPC). The spin system evo- 
lution will be an interplay of the non-unitary dynamics 
due to measurement and the dynamics associated with 
the system's internal Hamiltonian. 

The purpose of this work is to develop the theory of 
quantum evolution of a DQD two-electron spin system 
under the continuous measurement by a linear charge 
detector such as a QPC that is based on the Bayesian 
approach i 9 ' 18 ' 19 Initially motivated by experiment 1 ^ we 
pose the question: "How the measurement by a QPC 
will affect the spin qubit state and is it possible to use 
continuous measurement for manipulation of the quan- 
tum state ?" The system state is continuously updated 
according to a given measurement record due to quan- 
tum back action^ On typical system's times when many 
electrons pass through the QPC, the quantum evolution 
is conditioned on the fluctuating detector current J(t). 

We identify regimes when the system-detector coupling 
becomes large relative to a typical system frequency, that 
leads to suppression of coherent tunneling transitions. 
The associated stabilization of the states in continuous 
measurement (see, e.g., Ref j 18 i 20 ' 21 ) is an analog of the 
well known quantum Zeno effecti 22 i 23 ' 24 Since the system- 
detector coupling is finite, the Zeno stabilization does not 
last forever. The stabilized states switch to each other 
with a rate T sw -C T. The switching happens on a time 
scale 1/r. In fact, on this time scale the switching tran- 
sition becomes irreversible. For an ensemble averaged 
description of the state evolution we will show that the 
irreversibility corresponds to an exponential decay (com- 
pare with Ref.— ) with the decay rate T sw . 

Correspondingly, in the case of a spin qubit we show 
via numerical simulation of the measurement process that 



the switching rate T sw plays the role of the spontaneous 
decay rate T sp for the two-level atom described above. 
Despite the quantum back action noise reflecting the 
noisy output, the spin qubit state will evolve according 
to Eqs.(TT]), ([2]) in a spin-charge basis that corresponds to 
the Zeno stabilized states. In particular, in our analysis 
the negative-result evolution of the qubit subsystem is 
not postulated but emerges from continuous noisy mea- 
surement evolution and Hamiltonian evolution at a mi- 
croscopic level. 

The paper is organized as follows. In Sec. II we re- 
view/derive the Bayesian equations for weak continuous 
measurement of a DQD two-electron spin system. In Sec. 
Ill we investigate the measurement collapse scenarios for 
different coupling regimes and give qualitative and quan- 
titative explanations of the obtained results. In Sec. IV 
we estimate the detector-system coupling and other im- 
portant parameters that could be relevant for an exper- 
iment to confirm our results. We also comment on the 
role of various sources of decoherence. 



II. DQD TWO ELECTRON SYSTEM AND 
MEASUREMENT DEVICE 

A. Hamiltonian and time scales 

The system of two electrons confined in the DQD, and 
the quantum point contact as a charge sensing detector 
are shown in FigfT] External gate voltages form the po- 
tential profile of the DQD and are well controlled in the 
experiment! 1 - In particular, changing the energy differ- 
ence between the dots, <5, allows continuous tuning of the 
charge configuration between (1, 1) and (0, 2): when the 
electrons are in the right dot (0, 2) the ground state is 
the spin singlet 15(0,2)). The interdot tunneling cou- 
pling Q allows coherent transitions between the singlet 
spin-charge states |5(1, 1)) and |5(0, 2)) which are close 
in energy (for 6 = 0, the states |5(1,1)) and |5(0,2)) 
are in resonance). Transitions between the triplet states 
\T(1, 1)) and \T(0, 2)) are suppressed due to large energy 
mismatch: the highly energetic triplet states ^±^(0,2)) 
are well separated in energy from the singlet |5(0, 2)) 27 
(AEst ~ 400/ieV) due to tight confinement and on site 
exchange interaction. In the (1,1) configuration the spin 
qubit subspace is formed by the singlet |5(1, 1)) and the 
triplet state |Xb(l,l)) with zero spin projection on an 
external magnetic field B] the applied field removes the 
triplet degeneracy by splitting off the |T±(1,1)) states 
(splitting A± w 2.5/LteV for the GaAs DQD described in 
Refi 1 ^). Thus, the DQD spin system is considered below 
as a three state quantum system (qutrit). 

Introducing the short notations for the relevant three 
states: |1) = |5(0,2)), |2) = |5(1,1)), |3) = |T„(l,l)), we 
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write the system hamiltonia n 2 7 ' 28 : 

H DQD = -S\1){1\ + ^(\1){2\ + \2){1\) 

+es \2)(2\+s T ]3><3|. (3) 

The presence of the |T± i0 (0, 2)) states induces small ex- 
change energies £<j, Et which can be taken into account 
perturbativelyj^I in what follows we will neglect these en- 
ergies that are of the order of tt 2 / AEg? <§C ^- Other pos- 
sible terms in the Hamiltonian TLdqd (that are present 
for a general qutrit system 2 ^) are suppressed due to con- 
servation of spin. 

In the total Hamiltonian of system plus detector 

TL = TLdqd + Tidet + Tlint, (4) 
the (low transparency) QPC detector— is described with 

H det = E 4 c i + Yl E » c ^ + H( Tc « c ' + H - c -); (5) 

l u 1 >U 

here the operator cj (ct ) creates an electron in the lower 
(upper) lead of the detector (Fig[T]) and the tunneling T 
between leads is assumed energy independent. 

The system-detector interaction can be written in 
analogy with the qubit case; for a general qutrit mea- 
sured by a hypothetical detector that can distinguish all 
three states one formally writes the interaction Hamilto- 
nian as proportional to the two [SU (3)] generators, that 
are diagonal in the spin-charge basis: 

Wi„t=^[A(|l)(l|-|2)(2|)0 C t Ci 

l,r 

+B (|1) <1| + |2)(2| - 2 |3)(3|) ® eta + H.c] . (6) 

Introducing detector transmission probability Tk for each 
system state \k), the average currents can be expressed 
as Ik = 2TTTkpip u e 2 V/U, k = 1,2,3, where V is the de- 
tector bias voltage, and pi and p u are the densities of 
states in the lower and upper lead. In the case of QPC, 
which cannot distinguish between states |2) = 1*5(1, 1)) 
and 1 3) = |7o(l, 1)) the tunneling amplitudes A and B 
are not independent: T 2 = \T-A+B\ 2 = T 3 = \T-2B\ 2 , 
since the corresponding average currents are equal, I2 = 
Is = For the state |1) = 1*5(0, 2)) the average cur- 

rent I(o.2) 7^ ^(i,i) an d it can be distinguished from the 
(1, 1) states. 

In what follows we consider a situation when the in- 
ternal detector dynamics given by Eqs.©, © is much 
faster than the two-spin dynamics due to TLdqd- Also, 
it is assumed that the voltage applied to detector is 
large so that a typical detector decoherence time is much 
smaller than the typical electron tunneling times; thus, 
ft/eV <C e//(o,2), e /I{i,i) "C The first inequality im- 
plies that coherences between different electron passages 
in the QPC can be neglected. Thus, the QPC detector 
will behave essentially classically on the typical time scale 
of the DQD two-spin dynamics. 



B. Continuous quantum evolution according to the 

result 

The quantum state evolution of an open quantum 
system^ 3 -! has an analog in a classical state estimation 
procedure^ it takes into account the actual measure- 
ment record that is imperfectly correlated with the sys- 
tem state. In a situation when the QPC is a weakly re- 
sponding detector, |AJ| = |J( 0)2 ) - J(i,i)| < I(a,2),I(i,i), 
when every tunneling electron brings a little information 
about the system state, it is reasonable to condition the 
system evolution on the quasicontinuous noisy detector 
current I(t). 9 Typical measurement time (that reflects 
accumulation of a signal-to-noise ratio of order one) can 
be introduced^, r meas = 2S /(AI) 2 ~ 1/T, where 
So — 2e/(i.i) is the low-frequency spectral density of the 
detector shot noise. [We have neglected small differences 
between shot noises.] The finite measurement time im- 
plies that given the noisy current I(t), the system state 
is updated gradually. Below we sketch the derivation of 
the measurement evolution equations for a qutrit related 
to the "informational" Bayesian approachj2i2i 

For the measurement evolution alone, the most 
straightforward way is to consider the elementary act 
of scattering of an incoming electron, \i), off the QPC 
tunnel barrier that depends on the charge state of the 
system, \k). The scattered state is expressed as a lin- 
ear transform of the initial state: If) = Sk\i) where the 
scattering matrix is (see, e.g., Refl35l) 

and r/j, tk are the reflection and transmission ampli- 
tudes. Following Jordan and Korotkov^ every tunnel- 
ing electron can be mapped to an ancilla qubit, whose 
basis states are the scattering states: \R) for a reflected 
electron, and \T) for a transmitted electron. Scatter- 
ing will entangle the states of the ancilla qubit and the 
qutrit. Then projective measurement on the ancilla will 
lead to POVM (positive operator-valued measure) mea- 
surement operators^ of the form: Mr = diag{ri, r 2 , r 3 } 
and Mt = diag{ii, £3} that satisfy the completeness 
condition M R M R + M T = 1. The diagonal form of 
these operators in the qutrit basis (the spin-charge basis) 
follows from the diagonal form of the interaction Hamil- 
tonian Tli n t, Eq.©. For collecting electrons in the upper 
lead, the count of an electron updates the qutrit density 
matrix according to a POVM formula p — MtpM^/Vt 
where Vt — Tr[Mj, Mxp] — J2k PkkTk is the total proba- 
bility to find an electron in the upper lead and Tk = \tk\ 2 
are the transmission probabilities introduced in the pre- 
vious section. If an electron was not counted a similar 
update takes place using the measurement operator Mr 
and the reflection probabilities, 1 — Tk = \rk\ 2 - In the 
basis where Mr,x are diagonal these evolutions take the 
form of Bayesian updates^ 

The evolution rate of the system density matrix will 
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be related to the rate of tunneling through the QPC. In- 
troducing an average number A St of (independent) tun- 
neling attempts per time St, the average current^ (given 
the system state \k)) is then I k — eAT k . For times of the 
order of individual tunneling times, St ~ e/I k , we can 
consider a negative-result evolution of the qutrit density 
matrix. For the evolution of the diagonal elements in 
case of no tunneling we write: 

Pkk(St) = Pkk {0) (1 - T k ) ASt /Norm. (8) 

(Norm is a proper normalization.) For a low trans- 
parency QPC, Tfc <C 1, the number of attempts is large, 
ASt ~ 1/T fc > 1 so that (l~T k ) ASt sa e ~ T » ASt = e~^ St . 
Thus, we obtain a negative-result evolution similar to 
Eq.([T]) where the spontaneous decay rate r sp is replaced 
by the tunneling rates — , each for every k = 1,2, 3. The 
evolution of the non-diagonal elements can be shown to 
be of the form similar to Eq.([2]). 

The above evolution can be hardly observed, since for 
typical system times, I/O, lyT 3> e/Ik, many electrons 
pass through QPC. For N independent attempts one 
can derive the conditional probability for m successful 
tunnelings— (per time r) given the system is in state 

V(m, N\k) = C%TF(1 - T k ) N ~ m , (9) 

so that the Bayesian update of the system (qutrit) den- 
sity matrix can be shown to be of the form: p kk — 

PkkV(m,N\k)/Norm and p' M = Pki\j P kk p\J PkkPu- 

We mention that the update related to Eq.©, is just 
a simple composition of N elementary updates, each cor- 
responding to tunneling (1) or no tunneling (0) of a 
single electron: e.g., for a result Q = (0,1,..., 0) the 
measurement operator Mq is just a multiplication of 
N diagonal operators, Mq = Mr 1 Mt 2 ■ ■ ■ Mr n . Also 
note that Eq.([9]) sums up over Cjy identical possibil- 
ities since QPC cannot distinguish different sequences 
Q\,Qi, ■ ■ ■ with the same total charge tunneled to the 
right lead. By representing the (random) number m 
of tunneling electrons per time interval t through the 
QPC current, m = rl(t,r) — ff +T I(t')dt' , and using 
De Moivre-Laplace limit theorem (N ^> 1, T k fixed) we 
can replace V(m, N\k) by a Gaussian distribution for the 
measurement result I(t,r), given the state \k); P k (I) — 
V^T^cxp [-(I(t, t) - hfr/So], where I k ~ eMT k /r 
is the average current, and So is the shot noise spectral 
density for a low transparency and weakly responding 
QPC. The update of the system density matrix can be 
written again as a Bayesian (informational) evolution: 9 

Pki (t + r) = p k i (t) — j , (10) 

with the total probability of a particular result I given 
by P(T)=J2 k p kk (t)P k (T). 

Differentiating Eq. (fTU)l over r (at r — > 0) one obtains 
a stochastic equation for the qutrit state evolution (in 



Stratonovich forrn^ 9 -) . Taking into consideration the (in- 
ternal) Hamiltonian evolution and possible dephasing, 
one can write: 

+ w-/i)(/(t)-^)} 

- j\Hdqd, p]ki ~ lu Pki- (11) 

Here I(t) is the formal limit (at r — ► 0) of the observed 
detector signal I(t,r). For numerical simulations of a 
measurement one complements 9 - Eq. lTTj) by 

i{t) = Y J pkk{t)h+m (12) 

k 

that is consistent with the statistics of P(I); here £(t) is 
a white noise with a spectral density S$ = So. Eq. pip 
is of the same form as an analogous evolution equation 
for a system of N qubitsj^ here the summation is over 
the three qutrit states. The dephasing rates, 7m, are 
related to the detector non-idealityj2ii£ experimentally 42 
and theoretically 9 - ! 35 ' 40 a QPC is close to an ideal detec- 
tor. An SET is usually highly non- ideal 9 , however it may 
reach ideality close to 1 in the co-tunneling or Cooper 
pair tunneling regime.— 

It is worthwhile to note that the measurement evolu- 
tion and the Hamiltonian evolution enter in Eq. ljll]) in- 
dependently; they just reflect the measurement (POVM) 
and unitary postulates applied to the system at a coarse 
grained time t e/I k 3> Ti/eV where the evolution is 
noisy and quasicontinuous. In Sec. Ill we will show via 
numerical solutions of Eq. pip that the continuous mea- 
surement evolution and the Hamiltonian evolution inter- 
play non-trivially so that new effective (negative-result) 
measurement and Hamiltonian evolutions of the system 
arise at a larger time scale. 

C. Ensemble averaged evolution of the system 

A total ignorance of a particular measurement result 
I(t) corresponds to a situation when the QPC detector is 
considered just as a part of a (Markovian) environment 
surrounding the system. Correspondingly, the density 
matrix available to such an observer (denoted as (pm)) 
will be quite different from that described by Eq. pTj) . 
The density matrix (p(t)) can be related to p(t), Eq. fTTj) 
by a formal procedure of an ensemble averaging over pos- 
sible results I(t) at every time moment t, similar to the 
classical probabilities* 3 ^ The averaging can be performed 
(for a sufficiently small time interval r), e.g., by using the 
total probability of a particular result /, given by P(I) 
and then simply adding the Hamiltonian evolution. The 
result will be a standard master equation 44 

(pki) = -Ffc/ {Phi) - -tJHdqd, (p)}m (13) 
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with ensemble-averaged dephasing rates 

r kl = (i k -ii) 2 /4S + iki. 



(14) 



For a quantum limited (ideal) detector, j)-t = 0, the 
dephasing rates produced due to averaging are just 
(4--fi) 2 / 45 o which is the minimum allowed by quantum 
mechanics ,£ i 18 i 40 i 41 

The individual dephasings, 7^ 0, may be a conse- 
quence of partial ignorance of the measurement result 9 - 
and are parameterized by the detector ideality (efE- 
ciency) 77 (0 < 77 < 1): lkl = (tT 1 - 1) (h - Ii) 2 /4S , 
i.e., Tki = (Ik — h) 2 /45 , o?7. Other sources of decoherence 
of the DQD system will be discussed in Sec. IV. 



III. EMERGENCE OF NEGATIVE-RESULT 
EVOLUTION IN THE QUBIT SUBSPACE 

When the unitary evolution of the three-state system 
is taken into account one has to explore the full Bayesian 
Eq. lfTTjl . which, as a rule, does not provide simple so- 
lutions. In what follows, using Eq. ljlip . we will perform 
numerical simulation of the measurement process for var- 
ious regimes of the system-detector dynamics. The non- 
trivial interplay of the quantum dynamics can be seen if 
one compares the effects of measurement evolution alone 
vs. total evolution. The measurement alone tends to col- 
lapse the system (qutrit) in either the state |1) or to the 
qubit subspace {|2), |3)} leaving states |2) and |3) unre- 
solved. However, adding the continuous coherent mixing 
of |1) and 1 2) by Hdqd leads to effective resolution of the 
states 1 2) and |3), so that a continuous collapse happens 
either to the state |3) or to the remaining spin-singlet 
subspace. 

If the continuous collapse happens to the singlet sub- 
space, at small coupling T the system will perform quan- 
tum oscillations, that are weakly perturbed by the mea- 
surement. When the coupling becomes large the picture 
is qualitatively different: one reaches the regime of Zeno 
stabilization of the system's singlet states. The latter is 
characterized by long time intervals when the QPC cur- 
rent is either j.) or /(o,2)j interupted by rare switch- 
ing between them; the qutrit state is correspondingly 
1 2) = 15(1,1)) or |1) ee 15(0, 2)). If the system was ini- 
tially into the qubit subspace, the continuous collapse 
takes the form of a slow negative-result evolution of the 
qubit state until it eventually switches to 15(0, 2)). 

The negative-result evolution is a Bayesian evolution 
conditioned by the information that the "qubit did not 
switch" . This evolution emerges in the Zeno regime as a 
solution of the underlying Bayesian stochastic evolution, 
Eq. tfTTj) . It is important to note that the key feature 
for establishing a negative-result evolution is the irre- 
versibility of the switching event (compare with Reffl2h. 
Within the stochastic measurement evolution, the irre- 
versibility is consistent with an important property of 
Eq. (|ll[) that is seen in numerical simulations: given 
a measurement record I(t), two evolutions that start 



from different initial states, p^(0) or p( 2 )(0), will be- 
come undistinguishable 9 -^ at a time scale of the order or 
greater than 1/T. In other words the system "forgets" 
its initial state, p(0) so that further evolution is domi- 
nated by the result, I(t) itself. The regime of current 
stabilization, (that is essentially classical^) then would 
correspond to irreversibility of the switching event. 



A. Collapse scenarios 

For a numerical simulation of the measurement pro- 
cess we use Eq. fTTl) . supplemented with the relation 
for the current signal, Eq. (|12p . where we incorporated 
the current level degeneracy for the states \2) and |3) 
(^2 = h = I n order to understand the collapse 

scenarios it is instructive to look at the stochastic equa- 
tion for the relevant density matrix elements transformed 
from Eq. (|lip to their Ito form 3 - 9 - (below we use % = 1). 

pu = -nim/312+pn (1 - Pu) ww- £(*)) (15) 
2AI 

P33 = -P11P33 -5— £(*), (16) 
DO 

P22 = -pn - P33, and 
P12 = i S P12 + i — {pn - P22) - Ti2 P12 
- (2/>u - 1) pia ^ 

Do 

pi 3 = I d P13 - I — P23 - 1 13 Pl3 

-(2p n -l)p 13 ^£(i), 

DO 



n 2A1 ., , 

-l TT P13 ~ PU P23 £(*)■ 
I Do 



(17) 

(18) 
(19) 



Li 1) is the current difference between 



P23 



Here AI = I(o,2) — 
the charge subspaces that can be distinguished by the 
QPC. The system-detector coupling explicitly enters in 

Eqs.(TT7]),(rTHD as r 12 = r 13 = (A/) 2 /45 + 712 = r = 

(A/) 2 /45 r7 

In the absence of the Hamiltonian, TLdqd = 0, the 
equation for pu becomes a pure noise and decouples; it 
can be heuristically considered as a (position dependent) 
Wiener process for the restricted variable < p\\ < 1 
with a diffusion coefficient, ~ PnO- ~ Pll)i approach- 
ing a minimum (zero) at the endpoints, p\\ — or 
pn = 1. Just as in the single qubit case2^ it sug- 
gests that the endpoints are the two possible attrac- 
tors for a given realization of the measurement. In the 
qutrit case, the measurement also preserves the ratio 
P22I P33 = & and the measurement evolution happens 
on a ray ^33 = (1 — pn)/(l + a) in the physical tri- 
angle of the plane (pn,p 33 ) [0 < Pkk < l,k = 1,2,3]. 
Switching on the Hamiltonian mixes the density matrix 
components leaving 7)33 unaffected; geometrically this is 
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represented as a horizontal moving in the plane (pn, P33) 
that changes a [causing the system to change from one 
ray to another]. From this point of view it is easy to 
imagine how the new attractors become P33 = (a hori- 
zontal ray) or the point P33 = 1. Finally, it is clear from 
Eq. (fT6|) that if the system is in cither of the subspaces, 
{|1), 1 2)} or 1 3), it will remain there: neither Hamiltonian 
nor measurement evolution mixes the subspaces, which 
is true for any DQD parameters, S, f2. The decoupling 
of subspaces is a consequence of the spin blockade and 
indistinguishability of the states |2) and |3) by measure- 
ment. 

The decoupling implies that the system will continu- 
ously collapse to one of the subspaces under a weak con- 
tinuous measurement as seen in numerical simulations, 
Fig|2] Since the corresponding ensemble averaged equa- 
tion for P33(t) is simply d(p33(t)) / dt = (ensemble aver- 
aging implies just nullifying of the noise in Ito form of the 
equations^) , then the ensemble averaged density matrix 
element is conserved: (/^(i)) = /033(C)). Thus, it must 
be that the fraction of members of the ensemble that col- 
lapses to P33(t) = 1 (after a sufficient measurement time 
t) is just /333(C)), which means that the probability of col- 
lapse to either of the subspaces is ^33(0) and 1 — ^33(0) 
respectively (i.e., the usual probability rules apply). 
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FIG. 2: (Color online) Weak coupling regime, V <C SI, at zero 
detuning, (a) A realization of P33(t) evolution (mixed initial 
state) showing collapse to the spin-singlet subspace {|1), |2)} 
(dashed line); weakly perturbed coherent oscillations of pn(t) 
(thick line) into the subspace. (b) Qutrit purification due to 
measurement with an ideal detector, r\ — 1. 



B. System evolution at small coupling. Quantum 
oscillations and spin blockade 

In order to characterize the emergence of a negative 
result evolution from continuous noisy measurement we 
will first consider the case when the system may estab- 
lish quantum oscillations weakly perturbed by the mea- 
surement, i.e., when negative-result evolution is not yet 
established. In particular, this happens for small system- 
detector coupling, r <C Q while the quantum dots detun- 
ing S is also small. For simplicity, we first consider the 
case of an ideal detector, r\ = 1. 

An example of the three-state system evolution is 
shown on Figj2^. In this particular realization of the 
measurement process, while starting from a mixed ini- 
tial state (pn(0) = 0, p 22 (0) = p 33 (0) = 1/2, Pij (0) = 
0,i 7^ j) P 33{t) — > (dashed line on FigJ^), i.e., the 
system is continuously collapsed to the spin-singlet sub- 
space, {|1) = |5(0,2)),|2) = 15(1,1))} after a transi- 
tion time of the order of 1/r. On the same time scale, 
weakly perturbed quantum oscillations are established in 
this subspace [solid line on FigJ^Ja.) for pn(t)]. 

The oscillation scenario in case of collapse to the sin- 
glet subspace is easily understood. Once p3z(t) gets to 
zero, Eqs. ([T5)) ■ ([!?)) just become of the same form as that 
describing quantum oscillations for a single qubitJ^ For 
zero detunning, Repi2 decouples (Jiep^lt) — > 0), while 
pii(i) — P22{i) = cos[fii + <p(t)] and 2lmpi 2 (t) — sm[f^ + 
<p(t)] oscillate with a unity amplitude (for 77 = 1), with 
phase <p(t) that slowly diffuses in timeJ^ Correspond- 
ingly, the oscillating scenario can be distinguished by the 



average current (I{t))t = (I(o.2) +^(i,i))/2 or by the cur- 
rent spectral density Si(uj) = 2 J_ Kj(t) exp(iu>T)dT 

(with Ki(t) = (I(t)I(t + r)) t - (/(*))? being the cur- 
rent correlation function) . Since the evolution of P33 (t) is 
transitional in character (see Figj2j), it will not affect the 
long-time average Ki(t). Then using, e.g., the methods 
of Refl2ll. it can be shown that the detector power spec- 
trum will have the same form as in the one-qubit case^ 
for any DQD parameters, 5, Jl. In particular, in the weak 
coupling regime (for 5 = 0) the spectrum Si(u>) exhibits a 
Lorenzian peak at the Rabi frequency f2 with a signal-to- 
noise ratio of Atj, and a width T = (I(o,2) — I(i,i)) 2 /^SqT]. 

On (Fig[2b) the qutrit purity is plotted. [We have de- 
fined it as Pur(i) = (3Tr[p(i) 2 ] - l)/2, so that Pur = 1 
for a pure state and Pur = for a totally mixed state.] It 
is seen that the qutrit is eventually reaching a pure state 
[even though with a random phase, tp{t)\ for a time of 
the order of 1/r. The purification of the qutrit state is 
yet another demonstration of non-trivial interplay of dy- 
namics. Indeed, Hamiltonian evolution alone conserves 
purity, while a measurement alone leaves states |2) and 
1 3) unresolved. The purification of the state is due to an 
effective resolution of the states |2) and |3) and the fact 
that no information is lost with an ideal measurement. 

For a measurement with a non- ideal detector, rj < 1, 
simulations show (if collapse happened to the spin-singlet 
subspace) that the amplitude of quantum oscillations is 
less than 1 and fluctuates in time. Correspondingly, the 
average purity (Puv(t)) t saturates at some lower value. 
Using Eqs. p3)) - (fT^|) we have derived in the weak cou- 
pling regime (compare with Ref. 45 ) (Pur(t)) t — 1 + 
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(3/4)[l/2?7 - v/(l + l/2r/) 2 - 2] so that the qutrit state 
remains mixed; for small 77 the qutrit purity approaches 
1/4 + 3/277. 

If collapse happens to p33(t) — ► 1, the state will pu- 
rify even for measurement with a non-ideal detector. In 
this scenario the average current is I(x t i) and the power 
spectrum is flat, Si(uj) — Sq [the detector signal is just 

/(t) =/(!,!}+ £(*)]. 



C. Large coupling. Zeno stabilization and 
emergence of negative-result evolution 

We now turn to the case of relatively large system 
detector coupling, Y ^> ft; the DQD detuning energy, 
8, can be either small or large. Numerical simulations 
confirm that collapse scenarios remain the same in the 
large coupling regime consistent with our argumenta- 
tion in Sec. Ill A. However, in the strong coupling case 
the evolution qualitatively changes; the typical collapse 
time to either of the subspaces, {|1), |2)} or |3), becomes 
much longer than 1/r. Also, instead of quantum oscil- 
lations within the spin-singlet subspace {|1) = 15(0,2)), 
|2) = 15(1,1))}, associated with Q, we have relatively 
long stabilization of the system state (Fig. [3]) in one of 
the two states since the measurement is trying to localize 
and "freeze" the system in a definite charge state. This 
is a manifestation of the quantum Zeno effect ! 22 ' 23 in the 
case of continuous measurement the detector is always 
coupled to the system, so the approach to quantum Zeno 
regime corresponds to the limit of stronger and stronger 
coupling. For a finite coupling the long stabilization peri- 
ods will be interrupted by rare switching events between 
the subspaces (the switching transition time is of the or- 
der of 1/r; see below). 

Anticipating the scenario of a negative-result evolution 
in the qubit subspace outlined in the beginning of Sec. 
Ill, we turn to the calculation of the switching rate be- 
tween the two singlet states, which will set the time scale 
of non-unitary dynamics. From an ensemble averaged 
point of view the qubit "decay" from |5(1, 1)) to |5(0, 2)) 
(triplet state |T (1, 1)) is spin blocked) is described as a 
smooth decrease of (pii(t)) and therefore the switching 
rate would be most easily obtained if we consider the 
master equation (fT3"|) that follows from averaging of mea- 
surement evolution. 

Due to current level degeneracy, h = h = ^ nc 
dephasing rate r 2 3, Eg. ([14"]) . vanishes just because the 
detector cannot distinguish between the corresponding 
states (no information can be obtained and therefore no 
information can be lost). Introducing r 12 = T\% = Y = 
(A/) 2 /45o?7 the master equation reads 



(A 



■[Hdqd, (p)_ 



kl 



(20) 





tfl/2« 

FIG. 3: (Color online) Strong coupling regime, Y Q. Zeno 
stabilization in the qubit subspace, pn = 0, for a long time 
until it switches to pn = 1 (thick dashed line). The noisy 
qubit evolution of a mixed initial state shown for P22 — psa 
and |/923 1 vs. postulated negative-result evolution given by 
Eqs. (|23|l - (|25|) (long dashed and short dashed lines). 



In order to calculate the switching rate from the state 
|1) = |5(0, 2)) we start with the initial condition p\\ = 1. 
Using Eq. p0"|) one derives the small time evolution: 



1 



n 2 t 2 



(21) 



Notice, that there is no linear term in the expansion, 
Eq. pTj) , which tells us that the decay is not exponential 
at a small time scale. This fact makes quantum Zeno 
effect physics possible. The ^-coefficient turns out to be 
-(ft) 2 /4 = -Tr[p 2 Hl QD ]+Tr[jm D Q D pH DQD ], i.e. it is 
determined by coherent (Hamiltonian) evolution alone, 
consistent with the discussion in Refs. [22ll23ll24l . 

It is instructive to find an approximate solution for 
(pii(i)) considering in (f2TJj) terms proportional to as a 
small perturbation (fl <C Y, for strong coupling). To first 
order in the perturbation f2 we obtain for Im(/?i2(i)): 



Im(p 12 (t)) = e 



-rt 



[Repi2(0) sinSt + Imp 12(0) cos St] 



Y 2 + S 2 



2 +C 



-rt 



[8 sin St -Y cos St] L(22) 



Using the relation, (pu(t)) — pn(0) — Q J Q Im(p 12 (f )} dt' , 
we can obtain an approximate solution for (pu(t)) and 
its small t expansion coincides with first few terms 
of the expansion, Eq. (p?Tj) . However on a time scale 
t* > 1/r, as seen in Eq. ([22|) . the exponential terms 
drop out and one reaches an expansion that has a lin- 



ear term: (pn(t*)) ~ 1 



2(r 2 +<5 2 ) 



t* . Similarly, one 



can show that if one starts from an initial state in the 
{|2), |3)}-subspace then on a coarse grained time scale 



(022(f)) " P2 2 (0) 



1 



q r 
2(r 2 +^) 



t* 



(A3s(*)) is just con- 



served. By solving numerically the master equation (|20|) 
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for (pn), (/O22), Re(pi2), Im(pi 2 ) one confirms that the 
decay of the subspaces is indeed exponential for large 
times and the switching rate between the subspaces is 



r — 

- 1 - ,S7/! 



n 2 r 



2(r 2 + 5 2 ) 



(23) 



Note that the strong coupling limit when T O (<5 
is arbitrary) implies that T sw <C T, i.e. the subspace 
life time is much longer than both 1/T and l/tt imply- 
ing quantum Zeno stabilization.— The exponential decay 
at times t* > 1/r is a sign of the irreversibility of the 
measurement^ that appears as a switching event. 

Having calculated the switching rate between the spin- 
charge states we can write (postulate) an ansatz for the 
time evolution of p(t) according to a given result. Start- 
ing from the qubit subspace, for times t* ^ 1/r one will 
be able to discriminate between the two current values 
or /(o,2)) an d thus to distinguish whether the sys- 
tem has decayed (switched) to the third state or not. 
The conditional probability for the state |2) = |5(1, 1)) 
not to decay by time t is given by i"2(f) = exp(— T sw t), 
while analogous probability for |3) = |Tb(l, 1)) is P$(t) = 
1 (due to spin-blockade). Using the quantum Bayes 
rule, similar to the two- level atom, Eqs.UJ),([2]) [see also 
Eq. (fT0]) ], one can write the effective negative-result evo- 
lution of the spin qubit subsystem given that it did not 
decay by time t: 



022(f) 



P22 (0)^2(0 
Ptot(t) 



, /933(f) 



,033(0) 

Ptot(t) 



P23(f) = P23(0) 



' ggaft) 1033(f) 
P22(0) ,033(0) 



-i0(t,5,n,r) 



(24) 



(25) 



where the total probability not to decay is given by 
Ptotit) = P2 2 (0)P 2 (f) +P3 3 (0)P 3 (f), and 0(t,*,fi,r) is 
an accumulated phase (see below). 

Consistency of the informational (Bayesian) approach 
requires that the negative-result ansatz for the p(t)- 
evolution be reproduced by the underlying evolution, 
Eq. lfTTj) . according to the noisy record I(t). In Fig[3] 
we show the negative-result evolution for /022(f) — P33{t) 
and I /923(f) I defined by Eas.(|25 j) .l|51 jl .p5 }) . versus the den- 
sity matrix evolution generated through simulation of 
the noisy measurement process via Eqs. (ITT|) . (|12p . The 
(noisy) evolution in the qubit subspace is quite regu- 
lar, and the state eventually approaches |3) = |T (1, 1)) 
(if the system did not switch). One can see that the 
two evolutions well agree in the strong coupling regime; 
the agreement is already established at T > 5fL The 
agreement reveals a non-trivial property of the Bayesian 
stochastic evolution Eq. (fTTjl . It means that the postu- 
lated negative-result evolution given by Eqs. (f2"4")) . (f2"5)l . 
can be actually derived from Eq. (|lip in the Zeno regime, 
as an interplay of measurement evolution and Hamilto- 
nian evolution at an underlying microscopic level. 

The time scale of the negative-result evolution is set 
by the switching rate T sw , Eq. (l2"3"l) , in which the detec- 
tor rate T = (AI) 2 /4S v = (AI) 2 /4S + 712 is the total 



rate, including the measurement rate and the additional 
rate due to detector non-ideality. Via numerical simu- 
lations we have confirmed the dependence of T sw on r\. 
We note that while the rate 712 would lead to dephas- 
ing in the spin-singlet subspace {|1),|2)} [see Eq. ljPT)) ]. 
it affects the spin qubit coherently^ in the sense that 
Eq. (|2"5")) preserves the coherence ratio IP23I/VP22 Pzz- The 
reason is the indistinguishability of the states 15(1,1)) 
and |To(l;l)) by the measurement. Particularly, a pure 
state remains pure, which we have confirmed numerically. 
(A similar conclusion was drawn in Ref[T3 from a com- 
pletely different viewpoint.) The mixed state will gener- 
ally purify (see the discussion in Reffl9l). The final qubit 
(and qutrit) purification happens on the same time scale 
Tjto ^ r -1 as the collapse to the spin-singlet or spin- 
triplet subspaces. 

The update for the non-diagonal element P23 in ([25]) 
reflects not only the conservation of coherence, but in- 
cludes an accumulated phase <p(t, 5, Q, T), which remains 
undefined by the negative- result ansatz. Stochastic nu- 
merical simulations by Eqs. (fTTjl . (fT2")l show that the phase 
is linear in time even for small detuning, \5\ < CI: 

4>(t) = e^(S, ^ -T) t and vanishes for 5 — (while noisy, 
it stabilizes just for times > 1/r). The coefficient £33 
can be interpreted as an energy splitting between the 
spin qubit states |2) and |3) induced by the negative- 
result measurement in the presence of the localized sin- 
glet state, |1). For small detuning, 5 ~ £7, the energy 

splitting is small, |e|3 I "C lefiT*! — ifjrj ^ nas the same 
sign as the exchange splitting, , that would be in- 
duced perturbatively. 

The effective energy splitting can be derived if one 
compares ensemble averaging of the negative-result evo- 
lution given by Eqs. ([24|) . (|25p with the ensemble aver- 
aged evolution of the "original" stochastic equations (jlip . 
Since the negative-result evolution is the Bayesian evolu- 
tion at a coarse grained time scale t* > 1/r, its ensemble 
averaging must coincide with the evolution given by the 
master equation (|20p considered at times > 1/r. Aver- 
aging of the negative-result evolution is straightforward 
using the probability not to decay, Ptotit)- Starting from 
the qubit subspace it gives, e.g., for the non-diagonal ma- 
trix clement: 



(p 23 (f)> =/>2 3 (0)e- 



-r» ra t/2 



(26) 



On the other hand Eq. ([20]) gives for the ensemble aver- 
aged evolution of (/023(f)), (pis(f)): 



(/023(f)) =-^(/0i 3 (f)> 



(27) 



CI 



( Pl3 (t)) = 1 S ( Pl3 (t)) - 1 - (p 23 (t)) - r ( Pl8 (t)>. (28) 

It can be solved exactly, and for the initial values ,013(0) = 
0, /023(0), we obtain: 

(p 23 (f)) = P23(0)e" ft/2 ( cosh Clt + ^ sinh Clt J (29) 
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where f = T - i 8, fi = Vf 2 - S7 2 /2. Taking the strong 
coupling limit T 3> 57 ((5 is arbitrary), at times t* > 1/T 
when some contributions to Eq. (|29|) are exponentially 
suppressed [similar to Eq. l[22j) ]. we reproduce the evo- 
lution for (p23{t)}, Ea. tpp]) . with an energy splitting: 



1> = |S(0,2)> 



it approaches e v ^ 1 ~ jg for large 8. Eq. (|30)) is con- 
firmed with a very good accuracy by the results obtained 
through direct simulation of the stochastic Bayesian evo- 
lution equations (HI]). Eas.pg ]) .(|25 ]) .([5n ]) suggest that the 
spin-charge states under continuous strong measurement 
can be re-interpreted as the "new" energy states, if the 
quantum evolution is monitored at times t* > 1/f. 



n 2 s 



4(r 2 + s 2 ) 



r. 



'2T' 



(30) 



D. Zeno stabilization at small coupling and large 
detuning 

Interestingly, Zeno stabilization can take place even 
at small coupling, if the detuning is sufficiently large, 
r -C 57 <C 8. Qualitatively, this can be understood 
by kinematical reasons as illustrated on the spin-singlet 
Bloch sphere, FigdJ In the stochastic Bayesian equations 
(fTTj) . for large detuning 8 3> £7, the Hamiltonian evolu- 
tion is a fast rotation with £l sys = \/ 8 2 + 57 2 around an 
axis close to the z-axis. On the other hand, the mea- 
surement alone tries to localize the state to either of the 
poles (the singlet states) acting along the meridians. The 
angular velocity along the meridians is 57 <C 57 ss/s ; also 
the oscillation amplitude is ~ 5 -C 1, so that the Hamil- 
tonian evolution is effectively suppressed due to averag- 
ing. Thus, even though the coupling is small, r <§; f2, it 
might be large with respect to an effective evolution rate 
along the meridians which quantifies a quantum Zeno ef- 
fect. Numerical simulation of noisy measurement evolu- 
tion shows that the DQD two-electron system either col- 
lapses to the spin-singlet subspace (similar to that shown 
in Fig(3j), or ends up at the triplet state. In the first sce- 
nario, the system is stabilized for a relatively long time 
3> r _1 in one of the singlet states while performing 
rare switchings between them. 

The switching rate T sw is given by the same formula, 
Eq. (f!?3")) . since the 57-terms in the master equation are still 
small perturbations with respect to the remaining terms 
as in Eq. . Solution of the master equation (f!?0"| con- 
firms that at a time scale t* > 1/T the decay from the 
subspaces is exponential as was in the case of a strong 
coupling, with the same switching rate. Correspondingly, 
numerical simulation of noisy measurement evolution of 
the DQD two-electron system, as shown on FigEK,b, 
is in a good agreement with the negative-result evolu- 
tion described by Eqs.(|24] ) .([25 |) . Generally, the agree- 
ment is established already at 8 > 357. Here, the initial 
qubit state was chosen to be a pure state and it remains 
pure as well as the purity of the total qutrit state (not 




2>=|S(1,1)> 



FIG. 4: (Color online) Bloch sphere in the spin-singlet sub- 
space illustrates the Zeno stabilization in case of weak cou- 
pling and large detuning. The projection of the Hamilton 
evolution (rotation around an axis close to z-axis) along the 
meridians (measurement evolution) is effectively small. 
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FIG. 5: (Color online) Evolution in the regime of small cou- 
pling and large detuning, (a) A realization when the qubit did 
not switch: pn(t) remains zero always. Evolution of P22 — psz 
shows collapse to the state |3) = |2o(l, 1)). (b) Oscillation 
of coherences, Rep23 and Imp23 (thin dashed and solid lines), 
and their envelop \p%z\ (thick solid line) vs. negative-result 
evolution of \p2s\ (dashed line) according to Eqs. 



shown). The negative- result evolution in the spin qubit 
subspace takes place as long as the system did not switch 
to the third state. In Fig[5]a realization of the measure- 
ment process is shown when the system did not switch 
at all. Indeed, in a situation when only a cr z -evolution 
is present in addition to the negative-result evolution, 
Eqs. (pM)) . (f!?5)) , the probability not to decay is given by 
Ptotit) = P3 3 {0) + P22(0)e~ rsmt . So, it approaches ,033(0), 
remaining finite for large times. 

The induced energy splitting leads to a a z evolution for 
Hep23, Im/?23, whose amplitude envelope is provided by 
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the negative-result evolution, Fig[5]b. For large detuning 
and small coupling the energy splitting is mainly due to 
perturbative influence of the third state: eflj — £23*"* = 
n 2 /4S as it follows from Eq.(H|. 



IV. AVAILABLE EXPERIMENTAL 
PARAMETERS 

1. Coupling strength and switching rate. The pos- 
sibility to detect single electronic charges in quantum 
dots (QD) via quantum point contact or single elec- 
tron transistor was recently demonstrated in experi- 
ments of different groups i 15 i 16 ' 17 The detector coupling, 
r = (A/) 2 /ASqT], essentially determines the measure- 
ment time T meas ~ 1/r to reach a signal-to-noise ratio 
of order one, where AI is the current difference signal in 
the detector due to the presence (or absence) of an extra 
electron charge e in the QD. For the charge sensitivity of 
an rf-SET 5q = 2.4x 10~ 5 e/VWz, provided in RefH3, one 
estimates rf,^ T ss 10 7 s _1 . Recent experiment of Rim- 
berg groups demonstrated a charge sensitivity of an rf- 
SET of 6q = 2.4 x 10~ 6 e/VWz. Due to quadratic depen- 
dence on sensitivity, T cx (5q)~ 2 , this amounts to a two 
orders of magnitude improvement^ rf^, T ~ 10 9 s _1 . 
For an rf-QPC^ we estimated Tg p p c ~ 10 6 s -1 . The 
high detector ideality of a QPC also may be lost in the 
rf-regimei^ 5 - 

The above estimates show that, in principle, both a 
weak coupling as well as a strong coupling regime of a 
negative-result evolution are experimentally reachable. 
For a typical tunneling^ t c ~ 10 \xeV the characteristic 
frequency f2 = 2t c /h » 3 x 10 10 s _1 is large. Given the 
range of DQD detuning, \8\ ~ — 300 fieV, one can reach 
a quantum Zeno stabilization if a relatively large detun- 
ing is taken. Smaller values of tunneling, t c w 1 [ieV 
or even t c « 0.1 fieV, are also reachable 26 so that a 
strong coupling regime may be realized. Typical values 
of switching rate T sw for the presented parameters are in 
the range 2 x 10 5 - 5 x 10 7 s _1 . 

2. Decoherence due to charge fluctuations. The de- 
coherence in the DQD two-electron system has various 
sources. The coupling to uncontrollable detector de- 
grees of freedom leads to detector non-ideality, 77, par- 
tially discussed in Section II. Another mechanism of de- 
coherence is due to systems's coupling to background 
charge fluctuation s 27 ' 56 ' 57 ' 58 that will lead to a non- 
zero dephasing 723 in the spin-qubit subspace. It was 
argue d 27 ! 56 ' 57 that such dephasing will be well suppressed 
in the far-detuned regime 6 S> t c , where it may be 
723 ~ 10 3 — 10 5 s _1 . However, close to the charge degen- 
eracy, 5 < t c , the dephasing may strongly increas o 57 ' 58 to 
723 ~ 10 6 — 10 7 s _1 due to higher sensitivity to fluctua- 
tions of the DQD parameters. In the latter case it will be 
difficult to see the effect of Zeno stabilization since 723 
becomes comparable to the switching rate T sw . 

3. Decoherence due to phonons. Yet another decoher- 
ence mechanism (also related to charge degrees of free- 



dom) is due to coupling to a phonon environment i 59 i 60 
Physically, the relevant process is the double-dot in- 
elastic tunnelin g 15 ' 26 from state |2) = \S(1, 1)) to state 
|1) = 15(0, 2)) quantified by the inelastic rate Ti n (S). Re- 
laxation process associated with the inelastic tunneling 
may lead to a contribution to the switching rate T sw , 
assuming a weak environment, Ti n <C f2,T. Associated 
contribution to the dephasing 712 (also of the order of 
rj n ) is expected not to destroy the negative-result evolu- 
tion in this case. 

Estimations of the inelastic rate2Z give the range of 
0.01 — lOOneV (corresponding to approximately 10 4 — 
10 8 s _1 ) depending on the energy splitting AE between 
the relevant states. Due to generic cubic dependence on 
the splitting, Ti n cx (AE) 3 , one can hope to find a range 
of not too large detuning 5 so that Ti n is in the range 
of 10 4 — 10 5 s _1 . Thus, eventually one may expect to 
manage the inequality Ti n <C T SWl so that the Zeno sta- 
bilization will effectively "fight" against decoherence. We 
note that various models of boson environment may af- 
fect the predicted negative-result evolution in different 
ways that deserve a separate investigation. 

4. Implications of the QD nuclei. The contact hyper- 
fine interaction of the electron spins with the surrounding 
nuclei spins in the DQD leads to entanglement with the 
uncontrolable spin-bath degrees of freedom ; 61 ' 62 ' 63 ' 64 in 
quasi-classical language it can be described as an effect of 
"inhomogeneous broadening" quantified by the random 
nuclear Overhauser field. 27 ' 65 ' 66 ' 73 The dephasing caused 
by such effects in the spin qubit have been experimentally 
measured 1 "^ (dephasing time is of the order of ~ 10~ 8 s); 
such strong dephasing may conceal any interesting quan- 
tum evolution. For a non-interacting bath it was shown 
that this type of decoherence can be completely elimi- 
nated by using various spin-echo techniques^ 7 - Realisti- 
cally, for a weakly interacting bath, application of such 
techniques would reduce the decoherence by a few orders 
of mag nitude.^ 7 -^ For a GaAs DQD spin qubit system 
a "true" decoherence time of the order of ~ 10~ 6 s was 
measured experimentally using simple spin-echoiii 5 - 

This suggests that spin-echo techniques could be of 
use also to reveal the negative-result evolution we are 
discussing in this paper. We leave for a future project 
the investigation of various possibilities for application 
of spin-echo techniques in conjunction with a negative- 
result measurement evolution in the context of manip- 
ulation and/or preparation of the state of a DQD spin 
qubit. 



V. CONCLUSION 

We have shown that a negative-result evolution of a 
spin qubit can effectively emerge out of noisy measure- 
ment of a DQD spin system by a linear detector such 
as a quantum point contact. The evolution emerges as 
an interplay of measurement (non-unitary) dynamics and 
Hamiltonian dynamics of the three-level system, when 
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quantum Zeno stabilization of the spin qubit subspace 
takes place. 

Besides implications to the theory of quantum mea- 
surements, our results may be of practical use for manip- 
ulation of a DQD spin qubit (for papers discussing quan- 
tum measurements as an important resource see, e.g., 
RefsEHlliiillOlzllIzlIzi)- Recent experiments on a 
single phase qubit provide an interesting manipulation 
of its state via negative-result measurement i 1 1 ' 75 These 
advances support the hope that experimental implemen- 
tation of negative-result evolution is also possible in a 



DQD spin qubit system. 



ACKNOWLEDGMENTS 

The authors would like to thank A.N. Korotkov, L.P. 
Pryadko, and A.J. Rimberg for fruitful discussions and 
remarks. This work at Ames Laboratory was supported 
by the Department of Energy — Basic Energy Sciences 
under contract No. DE-AC02-07CH11358. 



1 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L.M.K. Vandersypen, Rev. Mod. Phys. 79, 1217 
(2007) 

2 M. Pioro-Ladriere, T. Obata, Y. Tokura, Y.-S. Shin, T. 
Kubo, K. Yoshida, T. Taniyama, and S. Tarucha, Nature 
Physics 4, 776 (2008). 

3 E. A. Laird, C. Barthel, E. I. Rashba, C. M. Marcus, M. P. 
Hanson, and A. C. Gossard, Phys. Rev. Lett. 99, 246601 
(2007) 

4 J. Berezovsky, M. H. Mikkelsen, N. G. Stoltz, L. A. Col- 
dren, and D. D. Awschalom, Science 320, 349 (2008). 

5 F. H. L. Koppens, K. C. Nowack, and L. M. K. Vander- 
sypen Phys. Rev. Lett. 100, 236802 (2008). 

6 J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 
68, 580 (1992). 

7 H. J. Carmichael, An Open System Approach to Quantum 
Optics (Springer, Berlin, 1993). 

8 H. M. Wiseman and G. J. Milburn, Phys. Rev. Lett. 70, 
548 (1993); Phys. Rev. A 49, 1350 (1994). 

9 A. N. Korotkov, Phys. Rev. B 63, 115403 (2001); Phys. 
Rev. B 60, 5737 (1999). 

10 P. S. Epstein, Am. J. Phys. 13, 127 (1945); M. Renninger, 
Zeitschr. fur Physik 158, 417 (1960); R. H. Dicke, Am. J. 
Phys. 49, 925 (1981); D. T. Pegg and P. L. Knight, Phys. 
Rev. A 37, 4303 (1988). 

11 N. Katz, M. Ansmann, R. C. Bialczak, E. Lucero, R. Mc- 
Dermott, M. Neeley, M. Steffen, E. M. Weig, A. N. Cle- 
land, J. M. Martinis, and A. N. Korotkov, Science, 312 
1498 (2006). 

12 L. P. Pryadko and A. N. Korotkov, Phys. Rev. B 76, 
100503(R) (2007). 

13 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

14 J. Elzerman, R. Hanson, L. H. Willems van Beveren, B. 
Witkamp, L. M. K. Vandersypen, and L. P. Kouwenhoven, 
Nature 430, 431 (2004). 

15 J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 
Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and 
A. C. Gossard, Science 309, 2180 (2005). 

16 R. Hanson, L. H. Willems van Beveren, I. T. Vink, J. 
M. Elzerman, W. J. M. Naber, F. H. L. Koppens, L. P. 
Kouwenhoven, and L. M. K. Vandersypen, Phys. Rev. Lett. 
94, 196802 (2005). 

17 W. Lu, Z. Ji, L. Pfeifer, K. W. West, and A. J. Rimberg, 
Nature 423, 422 (2003). 

18 A. N. Korotkov and D. V. Averin, Phys. Rev. B 64, 165310 
(2001); A. N. Korotkov, Phys. Rev. B 63, 085312 (2001). 



19 R. Ruskov, A. Mizel, and A. N. Korotkov, Phys. Rev. B 
75, 220501(R)(2007). 

20 R. Ruskov and A. N. Korotkov, Phys. Rev. B 67, 241305 
(2003). 

21 R. Ruskov, A. N. Korotkov, and A. Mizel, Phys. Rev. B 
73, 085317 (2006). 

22 L.S. Khalfin, JETP Lett. 8, 63 (1968). 

23 B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 756 
(1977); C.B Chiu, E. C. G. Sudarshan, and B. Misra, Phys. 
Rev. D 16, 520 (1977). 

24 A. Peres, Am. J. Phys. 48, 931 (1980). 

25 A. N. Korotkov and A.N. Jordan, Phys. Rev. Lett. 97, 
166805 (2006). 

26 F. H. L. Koppens, J. A. Folk, J. M. Elzerman, R. Hanson, 
L. H. Willems van Beveren, I. T. Vink, H. P. Tanitz, W. 
Wegscheider, L. P. Kouwenhoven, and L. M. K. Vander- 
sypen, Science 309, 1346 (2005). 

27 J. M. Taylor, J. R. Petta, A. C. Johnson, A. Yacoby, C. 
M. Marcus, and M. D. Lukin, Phys. Rev. B 76, 035315 
(2007). 

28 O. N. Jouravlev and Yu. Nazarov, Phys. Rev. Lett. 96, 
176804 (2006). 

29 C. M. Caves and G. J. Milburn, Optics Communications 
179, 439 (2000). 

30 S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 

31 E. B. Davies, Quantum Theory of Open Systems, (Aca- 
demic, London, 1976). 

32 A. S. Holevo, Statistical Structure of Quantum Theory, 
(Springer, Berlin, 2001). 

33 A. Shnirman and G. Schon, Phys. Rev. B 57, 15400 (1998). 

34 A. N. Korotkov, Phys. Rev. A 65, 052304 (2002). 

35 D. V. Averin and E. V. Sukhorukov, Phys. Rev. Lett. 95, 
126803 (2005). 

36 A. N. Jordan and A. N. Korotkov, Phys. Rev. B 74, 085307 
(2006). 

37 M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information, (Cambridge University Press, 
Cambridge, 2000). 

38 Independent tunneling events amount to the requirement 
of high detector voltage as mentioned in Sec. Ill Al Number 
of classical tunneling attempts per second can be evaluated 
as A = fp^r, consistent with the Landauer formula for 
average currents, given just after Eq.©. 

39 B.Oksendal, Stochastic differential equations (Springer, 
Berlin, 1998). 

40 D. V. Averin, Fortschr. Phys. 48, 1055 (2000); H. S. Goan 
and G. J. Milburn, Phys. Rev. B 64, 235307 (2001); A. A. 



12 



Clerk, S. M. Girvin, and A. D. Stone, Phys. Rev. B 67, 
165324 (2003); A. N. Jordan and M. Buttiker, Phys. Rev. 
Lett. 95, 220401 (2005). 

41 R. Ruskov, A. N. Korotkov, and A. Mizel, Phys. Rev. Lett. 
96, 200404 (2006). 

42 E. Buks, R. Schuster, M. Heiblum, D. Mahalu, and 
V. Umansky, Nature 391, 871 (1998). 

43 A. B. Zorin, Phys. Rev. L ett. 76, 4408 (1996); D. V. 
Averin, |cond-mat/0010052| (unpublished); A. A. Clerk, S. 
M. Girvin, A. K. Nguyen, and A. D. Stone, Phys. Rev. 
Lett. 89, 176804 (2002). 

44 A. O. Caldeira and A. J. Leggett, Ann. Phys. (N.Y.) 149, 
374 (1983); W. H. Zurek, Physics Today 44, 36 (1991). 

45 Qin Zhang, R. Ruskov, and A. N. Korotkov, Phys. Rev. B 
72, 245322 (2005). 

46 From an ensemble averaged point of view, when the result 
of the QPC is ignored, any qutrit state will be "frozen" for 
times r _1 < t < r s „. 

47 A. Peres, Quantum Theory: Concepts and Methods 
(Kluwer Academic Publishers, The Netherlands, 1995). 

48 Generally, we expect that possible presence of dephasing 
in Eq. (|25p . e.g. a factor e ' off *, will be strongly suppressed 
in the Zeno regime: j c g ~ P 2 ™/r <C T sw (compare with 
Ref42). 

49 A. N. Korotkov, Phys. Rev. B 78, 174512 (2008). 

50 W. W. Xue, B. Davis, F. Pan, J. Stettenheim, T. J. 
Gilheart, A. J. Rimberg, and Z. Ji, Appl. Phys. Lett. 91, 
093511 (2007). 

51 There is a room for a factor of 2 — 10 improvement of 
the charge sensitivity 8q (depending on the SET opera- 
tional temperature^) until reaching the theoretical limit 
obtained in Refl53l 

52 H. Brenning, S. Kafanov, T. Duty, S. Kubatkin, and P. 
Delsing, J. of Appl. Phys. 100, 114321 (2006). 

53 A. N. Korotkov and M. A. Paalanen, Appl. Phys. Lett. 74, 
4052 (1999). 

54 D. J. Reilly, C. M. Marcus, M. P. Hanson, and A. C. Gos- 
sard, Appl. Phys. Lett. 91, 162101 (2007); M. Thalakulam, 
W. W. Xue, F. Pan, Z. Ji, J. Stettenheim, Loren Pfeiffer, 
K.W. West, and A. J. Rimberg, "Shot-Noise-Limited Op- 
eration of a Fast Quantum-Point-Contact Charge Sensor" 
cond-mat . mes- hall/0708 .0861 ( 2007) . 

55 N. P. Oxtoby, J. Gambetta, and H.M. Wiseman, Phys. 
Rev. B 77, 125304 (2008). 

56 W. A. Coish and D. Loss, Phys. Rev. B 72, 125337 (2005). 

57 X. Hu and S. Das Sarma, Phys. Rev. Lett. 96, 100501 



(2006). 

58 A. Romito and Y. Gefen, Phys. Rev. B 76, 195318 (2007). 

59 A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. 
Fisher, A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 
1 (1987). 

60 T. Brandes and B. Kramer, Phys. Rev. Lett. 83, 3021 
(1999); T. Brandes and T. Vorrath, Phys. Rev. B 66, 
075341 (2002). 

61 V. V. Dobrovitski et al., |quant-ph/0112053| W. Zhang, V. 
V. Dobrovitski, K.A. Al-Hassanieh, E. Dagotto, and B.N. 
Harmon, Phys. Rev. B 74, 205313 (2006). 

62 V. V. Dobrovitski, H. A. De Raedt, M. I. Katsnelson, B. 
N. Harmon, Phys. Rev. Lett. 90, 210401 (2003). 

63 F. M. Cucchietti, J. P. Paz, and W. H. Zurek, Phys. Rev. 
A 72, 052113 (2005). 

64 A. V. Khaetskii, D. Loss, and L. Glazman, Phys. Rev. Lett. 
88, 186802 (2002). 

65 I. A. Merkulov, A. L. Efros, and J. Rosen, Phys. Rev. B 
65, 205309 (2002); S. I. Erlingsson and Yu. V. Nazarov, 
Phys. Rev. B 66, 155327 (2002). 

66 V. V. Dobrovitski, A. E. Feiguin, D. D. Awschalom, and 
R. Hanson, Phys. Rev. B 77, 245212 (2008). 

67 W. M. Witzel, R. de Sousa, and S. Das Sarma, Phys. Rev. 
B 72, 161306(R) (2005); W. M. Witzel, and S. Das Sarma, 
ibid 74, 035322 (2006); W. Yao, R.B. Liu and L.J. Sham, 
ibid 74, 195301 (2006); N. Shenvi, R. de Sousa, and K. B. 
Whaley, ibid 71, 224411 (2005). 

68 W. Zhang, V. V. Dobrovitski, Lea F. Santos, Lorenza Vi- 
ola, and B.N. Harmon, Phys. Rev. B 75, 201302 (2007). 

69 W. Mao, D.V. Averin, R. Ruskov, and A. N. Korotkov, 
Phys. Rev. Lett. 93, 056803 (2004). 

70 C. W. J. Beenakker, D.P. DiVincenzo, C. Emary, and M. 
Kindermann, Phys. Rev. Lett. 93, 020501 (2004). 

71 H.-A. Engel and D. Loss, Science 309, 588 (2005). 

72 A.N. Jordan, B. Trauzettel, G. Burkard, Phys. Rev. B 76, 
155324 (2007). 

73 W. X. Zhang, N. Konstantinidis, K. A. Al-Hassanieh, V. 
V. Dobrovitski, J. of Phys. - Cond. Mat. 19, 083202 (2007). 

74 X.-B. Wang, J.Q. You, and F. Nori, Phys. Rev. A 77, 
062339 (2008). 

75 N. Katz, M. Neeley, M. Ansmann, R. C. Bialczak, M. 
Hofheinz, E. Lucero, A. O'Connell, H. Wang, A. N. Cle- 
land, J. M. Martinis, and A. N. Korotkov, Phys. Rev. Lett. 
101, 200401 (2008). 



